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The equation of state and the phase diagram in two-flavor QCD are investigated by the Polyakov-loop ex- 
tended Nambu-Jona-Lasinio (PNJL) model with an entanglement vertex between the chiral condensate and 
the Polyakov-loop. The entanglement-PNJL (EPNJL) model reproduces LQCD data at /i > better than the 
PNJL model. Hadronic degrees of freedom are taken into account by the free-hadron-gas (FHG) model with 
the volume-exclusion effect due to the hadron generation. The EPNJL-l-FHG model improves agreement of the 
EPNJL model with LQCD data particularly at small temperature. The quarkyonic phase survives, even if the 
correlation between the chiral condensate and the Polyakov loop is strong and hadron degrees of freedom are 
taken into account. However, the location of the quarkyonic phase is sensitive to the strength of the volume 
exclusion. 

PACS numbers: 11.30.Rd, 12.40.-y 



I. INTRODUCTION 

A key issue in the thertnodynamics of quantum chromo- 
dynamics (QCD) is whether the chiral-symmetry restoration 
and the confinement-to-deconfinement transition coincide or 
not. If the two transitions do not coincide, phases such as 
constituent quark phase [Tl, 13] or quarkyonic phase may 
appear. When the two transitions are first order, discontinu- 
ities can appear simultaneously in their (approximate) order 
parameters, chiral condensate a and Polyakov loop HIM]- 
At zero quark chemical potential (/i), however, the two tran- 
sitions are found to be crossover by lattice QCD (LQCD) ^- 
[Tlll . Hence, there is no a piioii reason why the two transitions 
coincide. Actually, there is a debate lfT2ll on the coincidence 
in LQCD data at zero ^ jgjjTlll. 

LQCD has the sign problem for finite fj, 81 3il . This is an 
important problem to be solved in future. Fortunately, LQCD 
data are available at imaginary /i 114-20^ and real and imag- 
inary isospin chemical potential [21-231 because of no sign 
problem in the regions. The data show that chiral and decon- 
finement crossovers coincide within the numerical accuracy. 
This coincidence indicates that there exists a strong correla- 
tion (entanglement) between a and (p. 

An approach complementary to first-principle LQCD with 
the sign problem is to build an effective model consistent with 
LQCD data and apply the model to the real-/i region. The 
Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model 
is designed to treat both the chiral-symmetry restoration and 
the deconfinement transition 1124143811 . For the thermal sys- 
tem with temperature (T) and imaginary chemical potential 
/i = iOT, LQCD has the periodicity of 27r/3 in 0. For temper- 
atures higher than some critical temperature Te, LQCD has 
the first-order phase transition at 6* = 7r/3 mod 27r/3. These 
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periodicity and transition were first proposed by Roberge and 
Weiss 13911 , and then called the RW periodicity and the RW 
phase transition, respectively. The PNJL model can repro- 
duce LQCD data for the RW periodicity and the RW phase 
transition ll5ol - [35ll . but not for the coincidence between chi- 
ral and deconfinement crossovers fi^. A current topic at 
imaginary /i is the order of the RW transition at the endpoint 
T — Te- Recent LQCD shows that the order is first-order for 
small and large quark masses, but the order is weakened and 
could be second order at intermediate masses lfT6l [itIi . The 
PNJL model cannot reproduce the property ll40t . 

In order to solve these problems in the PNJL model, we re- 
cently extended the model so as to have a four-quark vertex 
depending on ^. If the gluon field has a vacuum expec- 
tation value (Aq) in its time component, such a four-quark 
vertex depending on # is allowed through (Aq) ||4i1i . The ver- 
tex introduced newly induces a strong correlation (entangle- 
ment) between <t and <P. The strong correlation is supported 
by recent calculations ll42[l of the exact renormalization-group 
equation |431. The new model is called the entanglement- 
PNJL (EPNJL) model |40]. The EPNJL results are consistent 
with LQCD data at zero and imaginary /i and also at finite 
isospin-chemical potential ^40]. Particularly, the model can 
reproduce not only the coincidence between chiral and decon- 
finement crossovers without ll40li and with the strong magnet 
field 1 44] but also the quark-mass dependence of the order of 
the RW endpoint mentioned above ll40ll . 

Another cuiTent topic related to the correlation between 
the chiral and deconfinement transitions is the quarkyonic 
phase ]3-6]. The concept of the phase is originally con- 
structed in the limit of a large number of colors (iVc). It is a 
confined phase with finite quark-number density (riq). Since 
the nq-generation is strongly induced by the chiral restoration, 
the quarkyonic phase is nearly equal to the chirally-symmetric 
confined phase. Recently, the PNJL calculations showed that 
the quarkyonic phase survives even for small Nc such as 3 IH- 
1^. It is then interesting whether the phase can survive in the 
EPNJL model with the strong correlation between a and <P. 
Another point to be argued is the relation between the quarky- 
onic phase and hadron degrees of freedom, since the location 
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of the phase was always discussed without thinking the de- 
grees of freedom for the case of — 3. 

In this paper, we analyze the thermodynamics of two-flavor 
QCD with the EPNJL model and its simple extension. First, 
we test the reliability of the EPNJL model by comparing the 
model results with LQCD data on the equation of state (EOS) 
at /i > 0. After confirming the reliability of the EPNJL model, 
we predict the phase diagram in the /i-T plane and investigate 
whether the quarkyonic phase survives under the strong cor- 
relation between a and Hadronic degrees of freedom are 
introduced by the free-hadron-gas (FHG) model with constant 
or T- and /i-dependent hadron masses. Here, the volume ex- 
clusion effect due to the hadron generation ||45[ l46ll is also 
taken into account in the FHG model. The volume exclu- 
sion effect guarantees the quark dominance in EOS at high 
T and/or high /i. EOS, the phase diagram and the presence or 
absence of the quarkyonic phase are investigated also by the 
EPNJLh-FHG model. 

In Section ini the PNJL model is recapitulated. In Section 
Hm the FHG model is explained briefly and numerical results 
of the EPNJL and EPNJLh-FHG models are shown. Section 
llVl is devoted to summary. 



II. PNJL MODEL 

We start with the standard two-flavor PNJL Lagrangian 
' 26 



C ^qi^i^uD" - mo)q 

+ G,[{qq)^ + {qtj^fqf] ~ U ^A] , <P[A]* , T) , (1) 

where q denotes the two-flavor quark field, toq represents the 
current quark mass, and D" = d'^ + lA'^ — iiiSq. Field A" 
is defined by A" = S^gA^- with gauge fields the Gell- 
Mann matrix A^, and the gauge coupling g. In the NJL sec- 
tor, Gs denotes the coupling constant of the scalar-type four- 
quark interaction. The Polyakov potential U, defined in (|7|, 
is a function of Polyakov loop <1> and its Hermitian conjugate 
(l>* 



1 



-TrL. <P* 



1 



TrL^ 



with 



L(x) —V exp 



(ir v44(x, r) 



(2) 



(3) 



where V is the path ordering and A4 — lA^. In the chi- 
ral limit (mo — 0), the Lagrangian density has the exact 
SU{Nf)i^ X SU{Nf)^ X f/(l)v X S'[/(3)c symmetry. The 
temporal component of the gauge field is diagonal in flavor 
space, because color and flavor spaces are completely sepa- 
rated in the present case. In the Polyak ov g auge, L can be 
written in a diagonal form in color space ll25ll : 



L : 



diag(e' 



',e 



3^/30c)^ (4) 



where = 03 + 08/\/3, 4>b = -(f'3 + (f>s/V^ and (j>c = 
— {(pa + Pb) = —2(j)s/V3- The Polyakov loop <P is an exact 
order parameter of spontaneous Z3 symmetry breaking in pure 
gauge theory. Although Z3 symmetry is not an exact one in 
the system with dynamical quarks, it still seems to be a good 
indicator of the deconfinement phase transition. Therefore, 
we use <!> to define the deconfinement phase transition. 

Making the mean field approximation and performing the 
path integral over the quark field, one can obtain the thermo- 
dynamic potential SI (per volume). 



i7 



-2Nf 



(27r)3 



3E{p) + i In J'gJ'g 



Um+U, 



Fq = 1 + 3(^ + <Z>*e-'^^ )e 

l + 3(r +<Z>e-^^+)e-^^^ 



,-3/3 -B^ 



(5) 
(6) 



where a 



(qq), M = mo - 2G's(t, Um 



Gsa^, E 



+ 1/2 and = E ± ^. On the right-hand side of dD, 
only the first term diverges. It is then regularized by the three- 
dimensional momentum cutoff A |4, 25, 26]. The variables 
X — <P,'P* and a satisfy the stationary conditions dSl/dX = 
0. 

We use U of Ref. HI, which is fitted to LQCD data in pure 
gauge flieory at finite T HS]: 



U = T^ 



a{T) 



+ b{T) ln(l - m<^* + 4(<?3 + <l>*^) - 3(^>^> 



(7) 



with 



(8) 



where the parameters are summarized in Table H] In pure 
gauge theory, the Polyakov potential yields a first-order de- 
confinement phase transition at T = Tq. The original value of 
To is 270 MeV that reproduces pure gauge LQCD data, but the 
PNJL model with the value of Tq yields a larger value of the 
pseudocritical temperature at zero chemical potential than 
Tc = 173 ± 8 MeV given by fufl LQCD data 1 49-51]. There- 
fore, we rescale To to 212 MeV in the PNJL model so as to re- 
produce LQCD result, T^. = 173 MeV, for the deconfinement 
transition temperature. However, the PNJL calculation yields 
the chiral transition temperature T„ — 216 MeV, while LQCD 
gives T„ = 173 MeV. Therefore, the PNJL model has a siz- 
able difference between T^r and T<|, , say Z\ = \Ta-T^\/T^ ^ 
20% m. Thus, the PNJL result is not consistent with LQCD 
data for and hence A. 

The sizable difference indicates that the entanglement be- 
tween the chiral and deconfinement transitions is weak in the 
PNJL model. In order to solve this problem, we proposed an 
effective coupling depending on the Polyakov loop, Gs(^). 
In fact, this vertex is discussed in the exact renormalization 
group method |41]. It is expected that ^ dependence of Gs (^) 
will be determined in future by an exact method such as the 
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OO 


ai 


£12 


63 


3.51 


-2.47 


15.2 


-1.75 



TABLE I: Summary of the parameter set in the Polyakov sector used 
in Ref. Ii27i1 . All parameters are dimensionless. 



exact renormalization group method 0411 - 14311 . In this paper, 
however, we simply assume the following G's(^) by respect- 
ing chiral symmetry, C symmetry ||34|| and extended Z3 sym- 
metry JH: 



(9) 



Thus, this model has entanglement interactions between a and 
<P in addition to the covariant derivative in the original PNJL 
model. The PNJL model with the entanglement vertex Gs(^) 
is referred to as entanglement-PNJL (EPNJL) model. The pa- 
rameters, ai and a2, determined from LQCD data at zero 
and imaginary chemical potentials are ai — a2 — 0.2 ll40ll . 
Hadron degrees of freedom will be considered in Section HID 



III. RESULTS 

A. Comparison of model calculation to LQCD data at fi > 

The equations of state (EOS) calculated with the PNJL 
and EPNJL models are compared with LQCD data at fi > 
ll52il53ll . Figure [T] shows the pressure p{T) and the energy 
density e{T) at /i = and the net quark-number density p{T) 



at fi ~ 0.8Tc, where p 



for the quark number 



density and the antiquark number density and is the 
pseudocritical temperature of the deconfinement transition at 
II — 0, that is, the peak position of the Polyakov-loop suscep- 
tibihty. Forp(T) and e(r), LQCD data [52J provide only the 
deviations, p(T) -p(T„) and e(T) -e{T^), from T„ = 0.9Tc. 
Hence, p{Tn) and e(rn) are evaluated by the free-gas model 
of hadrons with vacuum masses. This procedure is reliable 
at T = 0.9Tc, because p{T) and e{T) are dominated by the 
hadron components there; we will discuss this point later in 
subsection unci The p{T) and e{T) thus estimated from the 
LQCD data are shown by the dots in panels (a) and (b); note 
that in these panels the p{T) and e{T) are normalized by the 
values in the Stefan-Boltzmann limit at p = 0. In panel (c), 
p{T) is nondimensionalized by T^. 

For all of p, e and p, the EPNJL results (solid lines) 
are more consistent with LQCD data than the PNJL results 
(dashed lines), as shown in Fig. [T] The entanglement inter- 
action makes the chiral symmetry restoration faster, that is, 
it shifts Tcr down to T$ ||4C|1 . As a consequence of the shift- 
down property, the EPNJL model has rapid change in p{T), 
e{T) and p{T) with T and hence reproduces the sharp change 
of LQCD result better than the PNJL model. Both the EPNJL 
and PNJL models underestimate LQCD results at T < T^, be- 
cause these models have no hadron component in their p{T), 
e(T) and p{T). We will return this point later in subsec- 
tion HIICl 




Fig. 1: T dependence of (a) the pressure at /i = 0, (b) the energy 
density at /i = and (c) the net quark- number density at p/Tc = 
0.8. The pressure and the energy density are divided by the values 
in the Stefan-Boltzmann limit at p — 0, while the net quark-number 
density is divided by T'^. The solid (dashed) lines show the EPNJL 
(PNJL) result. LQCD data are taken from Ref. (52|| for panels (a) 
and (b) and Ref. L53.1 for panel (c). 



B. Phase diagram inT — p plane based on the EPNJL model 

The pseudocritical temperature of the deconfinement 
crossover was defined in two ways so far. One uses the peak 
of the Polyakov-loop susceptibility x<i> lUll^ and another does 
^ = 0.5 pS]. The two definitions yield almost the same value 
of T4, for lower p, but for higher p the former becomes some- 
what obscure since x<i> has a broad peak there. We then take 
the latter definition in this paper. 

Figure 12] (a) shows the phase diagram predicted by the EP- 
NJL model with constant Tq. The thin-solid line represents the 
chiral crossover defined by the peak of the chiral susceptibil- 
ity, while the dashed line shows the deconfinement crossover 
defined hy (p — 0.5. The thick-solid line stands for the first- 
order chiral transition. For p < 0.15 GeV, the entanglement 
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interaction makes the chiral and deconfinement crossovers al- 
most coincide. For /x > 0.15 GeV, however, the first-order 
chiral transition line and the deconfinement crossover line 
diverge, so that there appears a chirally-symmetric but con- 
fined phase between the two lines. Since the Tiq-generation is 
strongly induced by the chiral restoration, the chiral transition 
line corresponds to the rtq -generation line. In this sense, the 
chirally-symmetric confined phase is the quarkyonic phase in 
which quarks are confined but riq is finite. Thus, the entangle- 
ment interaction does not make the two transitions coincide 
for fj, > 0.15 GeV. This can be understood as below. 

When the entanglement interaction is switched off, in gen- 
eral, T<f differs from T^. When Tc, > T,p, is large at 
< T < Tcr, so that the four-quark interaction (|9) is sup- 
pressed there by the entanglement terms having ai and a2- 
This suppression shifts down to T<p. This is the situation 
for /i < 0.15 GeV. When < T^, meanwhile, (p is small at 
Ta <T < T,p, so that little entanglement effect occurs in (|9|l. 
This is the case for /x > 0.15 GeV. Thus, the entanglement 
effect takes place mainly when > T^. 
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Fig. 2: Phase diagram predicted by the EPNJL model with (a) con- 
stant To and (b) ^-dependent Tq. The thick (thin) solid line corre- 
sponds to the chiral transition line of first-order (crossover), while 
the dashed line represents the deconfinement crossover defined by 
$ — 0.5. The plus symbol stands for the critical endpoint of the first- 
order chiral transition. The hatching region labeled by 'Qy' stands 
for the quarkyonic phase of definition 1 . 

In principle, the Polyakov potential U may depend on /i as a 
consequence of the backreaction of the Fermion sector to the 
gluon sector In fact, /i dependence of Tq in U was estimated 
with the renormalization-group argument ||28| : 



for b{p) ^ 29/(67r) - 32n^/{7rT^) with ao = 0.304 and 
Tr = 1.770 [GeV]. Figure ^h) shows an effect of ^- 
dependent To on the phase diagram. The chiral transition 
line is not changed much by introducing the /^-dependent Tq. 
However, T4, lowers more for larger /i, since so does Tq. 
Eventually, the chirally-symmetric confined phase, i.e., the 
quarkyonic phase, between the first-order chiral transition line 
and the deconfinement crossover line shrinks. 

The concept of the quarkyonic phase |3] is originally con- 
structed in the large Nc limit. It is a confined (color-singlet) 
phase with finite ?iq. The definition of the quarkyonic phase 
becomes somewhat unclear for small iVc = 3, since the de- 
confinement transition is crossover there. The confined state 
is defined by ^ < 0.5, but the color-single state is by ^ = 0. 
Hence, a possible definition of the quarkyonic phase is a phase 
of riq 7^ and <P < 0.5. Another possible definition is a 
phase of nq 7^ and ^ that can exist only at T = 
and fi > M ^ 340 MeV The second definition seems to be 
too strict, when the deconfinement transition is crossover So 
we think the third possibility with the PNJL model. For this 
purpose, the net quark-number density p is divided into three 
components 



Pl+ P2+ P3, 



(11) 



where 



Pk 



2Nf 



(2^ 



k mfc(p) 



k mk(p) 



l + E,■"^J(p) l + E,-^i(p) 



(12) 



with 

mi = Me-^^-, 
7712 = 3^>e~2/3£;4 



mi = 3<?*e^'^^+ 



ma = 3<?*e"^ 
7713 = e~^'^^+. 



(10) 



Here, pk means the net quark-number density where k quarks 
or k antiquarks are statistically in the same state f^JS]. The 
first (second) term on the right-hand side of (ITZt represents 
the quark (antiquark) contribution. In order for a phase to be 
color-singlet, three quarks must have the same momentum. 
This indicates that the quarkyonic phase can be defined by the 
Pa -dominated region of p3 > pi,p2 flB- This statement 
is reliable as follows. For zero temperature, p is finite only 
when pL > M, that is, at p > 0.34 GeV. In the region of small 
T and p > 0.34 GeV, furthermore, p is dominated by pa be- 
cause of e^^ 3> 1 in (fTTT l. At p > M and low T where the 
Pa -dominated region emerges, the quark part is larger than the 
antiquark part in (fTZt . and the denominator 1 + J^j (p) °f 
the quark part is dominated by Za-invariant T^a- Thus, the pa- 
dominated region possesses the color-singlet nature approxi- 
mately. 

In Fig.[3j the quark phase appears out of the chiral transition 
line (the thin and thick solid lines). The quark phase is sep- 
arated into the pa- and pi -dominated regions by the dashed 
line, and the pa -dominated (pi -dominated) region is located 
below (above) the border, as expected. At p > M, the pa- 
dominated region is much wider than the chirally-symmetric 
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confinement region sliown in Fig.|2] If the quarkyonic phase 
is defined by the p3-dominated region, the phase appears at 
small T and ji > M, independently of 

Now we summarize two definitions of the quarkyonic 
phase. 

1. Phase of riq 7^ and <P < 0.5. 

2. Phase of p3 > pi,p2- 

The quarkyonic phase of definition 1 (2) is plotted by a region 
labeled by 'Qy' in Fig.|2](Fig.[3]i. For /i > M, the quarkyonic 
phase of definition 2 is wider than that of definition 1 . Under 
definition 2, the region of the quarkyonic phase is independent 
of ^, and also not so sensitive to the choice of constant Tq or 
/i-dependent Tq, as shown by the comparison of panels (a) and 
(b) in Fig. [3] 

Hereafter, we take definition 1 commonly used. We will 
show in subsections IIII CI and IIII Dl that the quarkyonic phase 
of definition 1 survives even after hadronic degrees of freedom 
are taken into account. Hence, so does the quarkyonic phase 
of definition 2. 




|i(GeV) 



Fig. 3: Phase diagram predicted by the EPNJL model with (a) con- 
stant To and (b) ^-dependent Tq. The thick (thin) solid line corre- 
sponds to the chiral transition line of first-order (crossover), while 
the dashed line represents a line of pi = ps. The hatching region 
labeled by 'Qy' stands for the quarkyonic phase of definition 2. 



C. Hadronic effect on phase diagram 

At small T and p, it is natural to think that QCD dynamics 
is governed by hadronic (mesonic and baryonic) modes than 
quark modes. A simplest way of treating these modes is the 



free-hadron-gas (FHG) approximation. Under this approxi- 
mation, the thermodynamic potential is just a sum of the PNJL 
thermodynamic potential i7qrk+giu, the free-meson thermody- 
namic potential i7insn and the free-baryon one i7bryn: 

m—7r,(T ^ ' 

(15) 

where = -^/p^ + and = -^p^ + ± 3/i for 
meson and baryon masses and mf,. Here, we consider 
proton and neutron as baryons and tt and a as mesons. 

It is well known that hadron has a finite volume. Once 
hadrons are generated in the thermodynamic system, some 
part of the volume of the system is occupied by generated 
hadrons. Owing to this volume-exclusion effect, the hadron 
generation is suppressed at high T and/or high p and hence 
the hadronic degrees of freedom do not contribute to EOS 
there. This is the property that EOS should satisfy at high T 
and/or high p, since hadrons are expected to disappear there. 
For this reason, the volume-exclusion effect has often been 
introduced to EOS. The volume-exclusion is described in the 
canonical ensemble with volume V and hadron number N by 
V — V — vN for the excluded volume v = |vrr'^ with ra- 
dius r. It is possible to convert the definition in the canonical 
ensemble to that in the grand-canonical ensemble by the fu- 
gacity and the Laplace transform ll45ll46ll . In this formulation 
of Refs. ||45l IMI, the total pressure P'^^'^^ with the volume 
exclusion effect is described by the total pressure P for point- 
particles as 

P--'\T,p,)=PiT,p,), (16) 

where the i stand for quark, antiquark, meson, baryon and 
antibaryon. In P, the chemical potential pi of ith particle 
species is replaced by the modified chemical potential 

p^=p^- V^P''^'' (17) 

with the excluded volume Vi — |7rrf of radius and {pi} — 
{p, —p, 0, 3/i, —3p) for quark, antiquark, meson, baryon and 
antibaryon, respectively. Thus, the volume-exclusion effect 
can be taken into account simply by replacing pi by pi in 
P. For simplicity, all hadrons are assumed to have the same 
volume. For r^, we take two cases; the charge radius of nu- 
cleon, 0.8 fm, and the radius of the repulsive core in the nu- 
clear force, 0.5 fm. The number density rii and the entropy 
density Si for the ith particle species with the volume exclu- 
sion effect are obtained from the thermodynamic consistency 
by 




(18) 
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with 



1 + 



(19) 



where hi and Sj are the number and the entropy density for 
point particles with the modified chemical potential ^i. This 
model is called the PNJL+FHG model in this paper In the 
PNJL+FHG model, the net quark-number density p is ob- 
tained by p = Uq — fiq + 3nb — Sfib for the baryon- and 
antibaryon-number densities rib and h\^. 

In this subsection, we assume that hadrons keep masses at 
vacuum even for the case of finite T and /i; namely, = 
139 MeV for pion, ma — 680 MeV for sigma meson and 
TOp = m,! = 940 MeV for proton and neutron. This ap- 
proximation is quantitatively valid at T < 0.9Tc, but even at 
T > 0.9Tc it is acceptable for qualitative discussion, as shown 
later in subsection llllDI 
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Fig. 4: T dependence of (a) the pressure at ^ = 0, (b) the energy 
density at /i = and (c) the net quark- number density at /i/Tc — 
0.8. The solid (dashed) lines represent results of the EPNJL+FHG 
(PNJL+FHG) model with the volume exclusion of r = 0.8 fm. See 
Fig.[T]for other information. 



Figure |4] represents the pressure, the energy density and 
the quark number density obtained by the PNJL+FHG and 
EPNJL+FHG models with constant hadron masses and the 
volume-exclusion of r = 0.8 fm. Comparing this figure with 
Fig. [H one can see for T < that the PNJL+FHG (EP- 
NJL+FHG) models give better agreement with LQCD data 
than the PNJL (EPNJL) models. For T > Tc, furthermore, 
the EPNJL+FHG model (solid line) yields a better coinci- 
dence with LQCD data than the PNJL+FHG result (dashed 
line). Thus, the EPNJL+FHG model well describes T depen- 
dence of LQCD data for all T up to 2Tc. 



> 
O 




o 



0.2 



0.1 



















M ' B \ 

1 


B 



0.2 0.4 

la(GeV) 



0.6 



Fig. 5: Dominant degree of freedom in fi-T plane. The quark- 
dominated region is out of the thick-solid line, while the meson- and 
baryon-dominated regions are located inside the thick-solid line. The 
meson- and baryon-dominated regions are separated by the thick- 
dashed line, and they are labeled by 'M' and 'B\ respectively. The 
hatching area labeled by 'Qy ' stands for the quarkyonic phase of def- 
inition 1. The calculation is done by the EPNJL+FHG model with 
(a) r = 0.8fm and (b) r = 0.5 fm. The thin-solid and thin-dashed 
lines correspond to the confinement and chiral transition lines, re- 
spectively. 

Figure |5] shows which mode is the main component of EOS 
in T — fi plane. The meson-, baryon- and quark-dominated 
regions are classified by comparing meson, baryon and quark 
degrees of freedom with each other; for example, the meson- 
dominated region is defined by the condition of Um > ?ib + 
and n,„ > riq + hq, where n,„ is the meson-number den- 
sity. The meson- and the baryon-dominated region are lo- 
cated inside the thick-solid line and they are separated by 
the thick-dashed line and labeled by 'M' and 'B', respec- 
tively. The results are obtained by the EPNJL+FHG model 
with r = 0.8 fm in panel (a) and r = 0.5 fm in panel (b). 
The meson-dominated region is located at small T and ^, and 
the baryon-dominated region is at low T and middle fi. The 
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bary on-dominated region is expanded by decreasing r, i.e., 
weaker volume-exclusion. 

In Fig. 121 the quark-dominated region is located out of the 
thick-solid line. The thin-solid and thin-dashed lines represent 
the confinement and chiral transition lines, respectively. The 
intersection between the quark-dominated region and the area 
between the thin-solid and thin-dashed lines is the quarkyonic 
phase of definition 1 ; this is shown by the hatching area with 
label 'Qy'. For stronger volume exclusion of r = 0.8 fm in 
panel (a), the quarkyonic phase is located at fi > 340 MeV 
and T < 100 MeV. For weaker volume exclusion of r = 
0.5 fm in panel (b), the quarkyonic phase shifts to higher /i, 
but an island of the phase remains around /i = 340 MeV and 
T = 120 MeV. 



D. Effect of T- and /i-dependent hadron mass on phase 
diagram 

In subsection nil CI we considered hadron degrees of free- 
dom by using the free gas approximation with constant mass. 
However, hadronic masses are changed with T and /i. In this 
subsection, the T and /i dependences are taken into account in 
a simple manner 

Meson modes are quantum fluctuations around the mean 
field and can then be calculated with the random phase ap- 
proximation (RPA). Up to order 1/iVc, the thermodynamic 
potential is obtained by llsil] 

f2 — i7MF + ^^RPA, (20) 

where J7mf is the mean-field part shown in (|5]l and J7rpa 
is the mesonic-fluctuation part described by the ring dia- 
gram iH: 

^2rpa = ? E / li^det[l - G.7J(g)] (21) 

with the mesonic polarization bubbles 

n,k{q) = -TY, J ^Tr[r;S{p + q)aS{p)] (22) 

for j, k = (J, 7r+, 7r_, ttq, where Tr is the trace in color, fla- 
vor and Dirac indices. S is the quark propagator in the mean- 
field approximation, S{p) — {'yi,{p — + M — /i7o) ^ ■ 
The meson vertex depends on meson taken; precisely, 
= 1, = iT+75, = iT_75, irs^s- 
Since it is difficult to calculate the dynamical mesonic fluc- 
tuations (l2Tl i exactly, we then make the pole approximation, 
that is, i^RPA is approximated into ]7,nsn of (fT4l i but with m„i 
replaced by the pole mass nij that satisfies 

det[l-Gsi7((?o ="b,q = 0)] =0- (23) 

Figure |6] shows T dependence of tt- and cr-meson masses 
by solid and dashed lines, respectively. The thick and thin 
lines correspond to results of the EPNJL and PNJL models, 
respectively. Since the chiral transition is sharp in the EPNJL 




0.8 1.2 1.6 



T/r, 

Fig. 6: Meson masses as a function of T at fi — 0. Solid (dashed) 
lines stand for pi (sigma) meson masses. The thick (thin) lines cor- 
respond to the EPNJL (PNJL) cases. 

model compared with in the PNJL model, T dependence of cr- 
meson mass is also changed rapidly in the EPNJL model. In 
the chiral-symmetry broken phase at small T and /i, vr-meson 
masses are small, so that i7 is dominated by J7nisn in its tt- 
meson part. In the chirally-symmetric phase at high T and/or 
high /i, on the contrary, tt- and cr-meson masses are getting 
large and hence the mesonic contribution to SI becomes small. 

The baryon mass is hard to obtain even with the RPA ap- 
proximation. In Ref. llssll . the baryon mass is calculated in the 
NJL model by regarding baryon as a bound state of a quark 
and a diquark. The baryon mass mi, is approximately de- 
scribed by the constituent quark mass M as 

nib w 2.24M + 0.18 GeV. (24) 

Here, we assume that J7biyn of ( fTSl l has baryon masses mt of 
(Ell. 

Now, we have SI with T- and /i-dependent rrim and mf,. 
Figure|7]is the same figure as Fig.|5] but T- and /i-dependence 
of hadron masses are taken into account in Fig. Q The vari- 
ation of hadron masses does not change qualitatively loca- 
tions of baryon-, meson-, quark-dominated regions and also 
a location of the quarkyonic phase; note that in panel (b) the 
quarkyonic phase emerges at fi much higher than 600 MeV. 
More precisely, the hadron-mass variation shifts the boundary 
between baryon- and quark-dominated regions to higher T, 
because of baryon-mass suppression near the boundary. 

IV. SUMMARY 

The equation of state in two-flavor QCD was investigated 
with the EPNJL model and its simple extension. The EPNJL 
result is consistent with LQCD data on EOS at /i > better 
than the PNJL model. Thus, the EPNJL model is more re- 
liable than the PNJL model. After confirming the reliability 
of the EPNJL model, we have predicted the phase diagram in 
the fi-T plane with the EPNJL model. The quarkyonic phase 
survives, even if the correlation between the chiral condensate 
and the Polyakov loop is strong. As an extension of the EP- 
NJL model, we have introduced hadronic degrees of freedom 
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0.2 0.4 
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by using the free-hadron-gas (FHG) model with constant or T- 
and /i-dependent hadron masses. Here, the volume exclusion 
effect due to the hadron generation is also taken into account 
in the FHG model. The volume exclusion guarantees that the 
quark degree of freedom is dominant in EOS at high T and/or 
high jjL. The EPNJL+FHG model improves agreement of the 
EPNJL model with LQCD data particularly at small T. The 
quarkyonic phase survives, even if hadron degrees of freedom 
are taken into account. However, the location of the quarky- 
onic phase is sensitive to the strength of the volume exclusion. 
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